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1. Introduction 

The unit 3-sphere centered on the origin is the set of IR 4 defined by 

§ 3 = {(xi,x 2 ,x 3 ,x 4 ) eR 4 : xj+xl + xj + xj^ 1)}. 

It is often convenient to regard R 4 as the two complex dimensional space C 2 or the space of quaternions 
EL The unit 3-sphere is then given by 

S 3 = {( Zl ,z 2 )eC 2 : \ Zl \ 2 + \z 2 \ 2 = I)} or S 3 = {q E H : \q\ 2 = 1)}. 

The latter description represents the sphere § 3 as a set of unit quaternions and it can be considered as 
a group Sp(l), where the group operation is just a multiplication of quaternions. The group Sp(l) is 
a three-dimensional Lie group, isomorphic to SU(2) by the isomorphism C 2 3 (2:1,2:2) — > q € H. The 
unitary group SU (2) is the group of matrices 

( Zl - I 2 ), zi,Z2€C, m 2 + m 2 = l 

\—Z2 Zi j 

where the group law is given by the multiplication of matrices. Let us identify M 3 with pure imaginary 
quaternions. The conjugation qhq of a pure imaginary quaternion ft by a unit quaternion q defines 
rotation in R 3 , and since \qhq\ — \h\, the map h qhq defines a two-to-one homomorphism Sp(l) — > 
50(3). The Hopf map tt : § 3 -> § 2 can be defined by 

S 3 3 q^qiq = ir(q) G § 2 . 

The Hopf map defines a principle circle bundle also known as the Hopf bundle. Topologically § 3 is a 
compact, simply-connected, 3-dimensional manifold without boundary. 
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Even a small part of properties of the unit 3-sphcrc finds numerous applications in complex ge- 
ometry, topology, group theory, mathematical physics and others fields of mathematics. In the present 
paper we give a new look at the unit 3-sphere, considering it as a sub-Riemannian manifold. The sub- 
Riemannian structure comes naturally from the non-commutative group structure of the sphere in a 
sense that two vector fields span the smoothly varying distribution of the tangent bundle and their 
commutator generates the missing direction. The sub-Riemannian metric is defined as the restriction of 
the euclidean inner product from R 4 to the distribution. The present paper devoted to the description 
of sub-Riemannian geodesies on the sphere. The sub-Riemannian geodesies are defined as a projection 
of the solution to the corresponding Hamiltonian system onto the manifold. We give explicit formu- 
las using different parametrizations and discuss the number of geodesies starting from the unity of the 
group. While working on this paper the authors became aware on the results in [3] where the Lagrangian 
approach was developed and the minimizers were found (in our terminology geodesies are solutions to 
a Hamiltonian system so a minimizer is one of them). 



2. Left-invariant vector fields and the horizontal distribution 

In order to calculate left-invariant vector fields we use the definition of § 3 as a set of unit quater- 
nions equipped with the following noncommutative multiplication "o": if x = (xx, Xi, x 3 , X4) and 
V = (2/1,2/2,2/3,2/4), then 

((xxyx - X2V2 - x 3 y 3 - x 4 y 4 ), 
{x 2 yi + xiy 2 - X4JJ3 + 2:32/4), 
(x 3 yi + x 4 y 2 + xiyz - a; 2 y 4 ), (2.1) 

(^42/1 



Xoy= (xi,X 2 ,x 3 ,x 4 ) o (2/1,1/2,2/3,2/4) 



£32/2 + £22/3 + a; l2/4 J 

The rule (|2.ip gives us the left translation L x (y) of an element y — (2/1,2/2,2/3,2/4) by the element 
x = (x%, X2, x 3 , X4). The left-invariant basis vector fields are defined as X(x) — (L x (y)) it X(0), where 
X(0) are basis vectors at the unity of the group. The matrix corresponding to the tangent map (L x (y))* 
calculated by (|2.ip becomes 

/ xi -x 2 —x 3 -x 4 \ 
X2 Xi —X4, x 3 
x 3 X4 xi -x 2 
\Xi -x 3 x 2 Xi ) 

Calculating the action of {L x (y))* in the basis of unit vectors of M 4 we get four vector fields 

Xi(x) = xid Xl + x 2 d X2 + x 3 d X3 + x 4 d X4 , 

X 2 (x) = -x 2 d Xl + xid X2 + x A d X3 - x 3 d Xi , (2.2) 

X 3 (x) = -x 3 d Xl - x 4 d X2 + xxd^ + x 2 d Xi , 

X 4 (x) = -x 4 d Xl + x 3 d X2 - x 2 d X3 + xxd Xi . 

It is easy to see that the vector Xi(x) is the unit normal to § 3 at x with respect to the usual inner 
product (•, •) in R 4 , hence, we denote Xi(x) by N. Moreover, 

(N,X 2 (x)) = (N,X 3 (x)) = (N,X±(x)) = 0, and \X k (x)\ 2 = (X k (x),X k (x)) = 1, 

for k = 2,3, 4, and for any x £ § 3 . The matrix 




has rank 3, and we conclude that the vector fields ^(x), X 3 (x), X^(x) form an orthonormal basis of 
the tangent space T X E> 3 with respect to (•, ■) at any point x £ S 3 . Let us denote the vector fields by 

X 3 =X, X A = Y, X 2 = Z. 
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The vector fields possess the following commutation relations 

[X, Y] = XY -YX = 2Z, [Z, X] = 2Y, [Y, Z] = 2X. 

Let V = span{X, Y} be the distribution generated by the vector fields X and Y, Since [X, Y] = 2Z ^ T>, 
it follows that T> is not involutive. The distribution T> will be called horizontal. Any curve on the 
sphere with the velocity vector contained in the distribution T> will be called a horizontal curve. Since 
T X E> 3 — span{X, Y, Z = 1/2 [X, Y]}, the distribution is bracket generating. We define the metric on the 
distribution T> as the restriction of the metric (•, •) onto T>, and the same notation (•, •) will be used. 
The manifold (§ 3 ,2?, (•, •)) is a step two sub-Riemannian manifold. 

Remark 1. Notice that the choice of the horizontal distribution is not unique. The relations [Z, X] = 2Y 
and [Y, Z] = 2X imply possible choices V = span{X, Z} or T> — span{F, Z}. The geometries defined 
by different horizontal distributions are cyclically symmetric, so we restrict our attention to the T> = 
span{X, Y}. 

We also can define the distribution as a kernel of the following one form 

u> = —X2olxi + x\dx2 + Xidxs — X30IX4 

on R 4 . One can easily check that 

LO(X) = 0, lu(Y) = 0, lu{Z) =1^0, lu(N) = 0. 

Hence, kerw = span{X, Y, N}, and the horizontal distribution can be written as 

S 3 3nP I = kerw n T X S> 3 . 

Let 7(s) = (xi(s), X2{s), X3(s), xa{s)) be a curve on S 3 . Then the velocity vector, written in the 
left-invariant basis, is 

7(s) - a(s)X( 1 (s)) + b(s)Y( 7 (s)) + c(s)Z( 7 (s)), 

where 

a = (7, X) = — xz±\ — X4X2 + xi± 3 + X2X4, 

b = (7, Y) = — X4X1 + X3X2 — X2X3 + xi±4, (2.3) 

C = (7, Y) = — X2X1 + Xi±2 + X4X3 — X3X4. 

The following proposition holds. 

Proposition 1. Letj(s) — (x\(s),X2(s),yi(s),y2(s)) be a curve onS 3 . The curve 7 is horizontal, if and 
only if, 

c = (7, Z) = (7, X) = -x 2 xi + xi±2 + X4X3 - X3X4 = 0. (2.4) 

The manifold § 3 is connected and it satisfies the bracket generating condition. By the Chow 
theorem j2], there exists piecewise C 1 horizontal curves connecting two arbitrary points on S 3 . In fact, 
smooth horizontal curves connecting two arbitrary points on S 3 were constructed in pQ. 

Proposition 2. The horizontality property is invariant under the left translation. 

Proof. It can be shown that (|2.3[) does not change under the left translation. This implies the conclusion 
of the proposition. □ 



3. Hamiltonian system 

Ones we have a system of curves, in our case the system of horizontal curves, we can define the length 
as in the Riemannian geometry. Let 7 : [0, 1] — > § 3 be a horizontal curve such that 7(0) = x, 7(1) = y, 
then the length £(7) of 7 is defined as the following 

1(7)= f\i,i) 1/2 dt = [ (a 2 (t)+b 2 (t)) 1/2 dt. (3.1) 
Jo Jo 
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Now we are able to define the distance between the points x and y by minimizing the integral (|3.ip or 
the corresponding energy integral J Q (a 2 (t) +6 2 (t)) dt under the non-holonomic constraint (|2.4p . This is 
a Lagrangian approach. The Lagrangian formalism was applied to study the sub-Ricmannian geometry 
of § 3 in [TJ [3] . In the Ricmannian geometry the minimizing curve locally coinsides with the geodesic, 
but it is not the case for the sub-Riemannian manifolds. Interesting examples and discussions can be 
found, for instance in [H [7J [H [5] ■ Given the sub-Riemannian metric we can form a Hamiltonian 
function defined on the cotangent bundle of § 3 . The geodesies in the sub-Riemannian manifolds are 
defined as a projection of the solution to the corresponding Hamiltonian system onto the manifold. It 
is a good generalization of the Riemannian case in the following sense. The Riemannian geodesies (that 
arc defined as curves with vanishing acceleration) can be lifted to the solutions of the Hamilton system 
on the cotangent bundle. 

In the present paper we are interested in the construction of sub-Ricmannian geodesies on (S 3 , T> , (■, •)). 
Let us write the left-invariant vector fields X, Y, Z, using the matrices 



h = 







1 
1 



h = 





1 

-1 





h = 



-1 

1 









1 

-1 



(3.2) 



H 



Then 

X = (hx,Wx), Y = {hx,Vx), Z = {I 3 x,Wx) 
The Hamiltonian function is defined as 

1 

2 l ~~ ' ' ' 2' 
where £ = Vx. Then the Hamiltonian system follows as 

. _ dH 

dH 

dx 

As it was mentioned, a geodesic is the projection of a solution to the Hamiltonian system onto the 
x-space. We obtain the following properties. 

1. Since (I\x,x) = (hx,x) = (hx,x) = 0, multiplying the first equation of (|3.3p by x we get 



x = (hx,£,) ■ (I lX ) + (I a x,Z) ■ {hx) 

t = (hx,o-(hO + (i 2 x,o-(hO. 



(3.3) 







const. 



We conclude that any solution to the Hamiltonian system belongs to the sphere. Taking the constant 
equal to 1 we get geodesies on § 3 . 
2. Multiplying the first equation of (|3.3[) by I^x, we get 



{x,hx) = 0, 



(3.4) 



by the rule of multiplication for I\ 1 and I3. The reader easily recognizes the horizontality 
condition (x,Z) — in (13. 4p . It means that any solution to the Hamiltonian system is a horizontal 
curve. 

3. Multiplying the first equation of (|3 . 3|) by I\x, and then by I2X, we get 

(€,hx) = (x,hx), {£,xl2) = (x,xl 2 ). 

From the other side, we know that (x,Iix) = a and (a;, xl 2 ) — b. The Hamiltonian function can 
be written in the form 

H = \ ((hx, 2 + (hx, 2 )=\ ((hx, xf + (I 2 x, x) 2 ) =\(a 2 + b 2 

Thus, the Hamiltonian function gives the kinetic energy H — ^4j- and it is a constant along the 
geodesies. 
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4. If we multiply the first equation of (|3 . 3[) by x, then we get 

\x\ 2 = (hx, 2 + {hx, 2 = (hx, x) 2 + (I 2 x : x) 2 =a 2 + b 2 = 2H. 

Therefore 



\x\ 2 = a 2 



b 2 . (3.5) 



4. Velocity vector with constant coordinates 

We know that the length of the velocity vector is constant along geodesies. Let us start from the simplest 
case, when the coordinates of the velocity vector are constant. Suppose that a = b = 0. The first line of 
system (|3.3j) can be written as 

x\ = — CLX3 — bxi X3 = +axi — bx2 (4-1) 

±2 = — ax 4 + 6x3 X4 =+ax2 + bx\. 

Differentiation of system (|4.1[) yields 

X\ = — a±3 — bx 4 £3 — +aii — bx 2 (4-2) 

X2 =—ax4 + bx3 £4 =+aX2 + bx\. 

We substitute the first derivatives from (|4.ip in (|4.2p . and get 

x k = -r 2 x kl r 2 = a 2 + b 2 , k = 1,2,3,4. (4.3) 

Theorem 1. The set of geodesies with constant velocity coordinates form a unit sphere § 2 in R 3 

Proof. We are looking for horizontal geodesies parametrized by the arc length and starting from the 
point x(0) = xo- So, we set r — 1 and a = cosip, b = sin ijj, where ^ is a constant from [0, 2ir). Solving 
the equation (|4.3p we get the general solution x(s) = Acoss + Bsins. We conclude that A — xq from the 
initial data. To find B let us substitute the general solution in equations (|4.1|) and get B = (ah+btyxo. 
Thus, the horizontal geodesies with constant horizontal coordinates are 

x(s) — xq coss + (cosi/jIi + sin?/'-^)^ sins 

Since the geodesies are invariant under the left translation it is sufficient to describe the situation at 
the unity element, e.g., x = (1, 0, 0, 0) of S 3 . In this case the geodesies are 

x\ =coss, X3 — cos ip sin s, (4-4) 

x 2 —0; X4 = sin tp sins. 

We see that the set of geodesies with constant velocity coordinates form the unit sphere S 2 in R 3 = 
{(x\, 0, £3, X4)}. The parameter ip € [0,27r) corresponds to the initial velocity. □ 

The sphere (|4.4p is a direct analogue of the horizontal plane in the Heisenberg group H 1 at the 
unity. We remark that this result was obtained independently in [3] . 

Let us calculate the analogue of the vertical axis in H 1 . We wish to find an integral curve for the 
vector field Z . In other words, we solve the system 

a = (j,X) — —X3X1 - x A x 2 + X1X3 + X2X4 = 0, 

b = (7, Y) = —X4X1 + X3X2 — X2X3 + X1X4 = 0, (4.5) 

C = (7, Z) = — X2X1 + Xi±2 + X4X3 — X3X4 = 1, 
n = { A f,N)=+XiXi+X2X2+X3X3+X4X4 = 0. 

The determinant of the system is 1 and it is redused to 

x\ — —x 2} x 3 = +x 4 , 

± 2 = +Xi, ±4 = —X 3 . 
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Differentiating again, we get the equation x — —x. The initial point is x(0) = xq. System (|4.5[) gives 
the value of the initial velocity i(0) = /3X0. Taking into account this initial data we get the equation 
of the vertical line as 

x{s) = Xq cos s + I3XQ sin s. 
In particular, at the point (1, 0, 0, 0) the equation of the vertical line is 

x-i = coss, X2 = sins, X3 = 0, X4 — 0, s £ [0, 2ir]. (4-6) 



5. Velocity vector with non-constant coordinates 

Cartesian coordinates 

Fix the initial point a;^°) = (1,0,0,0). It is convenient to introduce complex coordinates z = 
x\ + ix2, w = X3 + ixi, ip = £1 + and ip — £3 + 1^4. Hence, the Hamiltonian admits the form 
2H = \wip — ztpl 2 . The corresponding Hamiltonian system becomes 

z = w(wip — zip), z(0) = 1, 

w = —z(w(p — zip), w(0) = 0, 

(p — ^j{wLp — zip), (p{Q)—A — iB, 

ijj = —(p(w(p — zip), ip(0) = C~iD. 

Here the constants B,C, and D have the following dynamical meaning: w(Q) = C + iD, and B = 
— ?ii)(0)/2w(0). So C,D is the velocity and -B/V C 2 + D 2 is the curvature of a geodesic at the initial 
point. This complex Hamiltonian system has the first integrals 

zip — wip — C + iD, 
zip + wip — A — iB, 

and we have \z\ 2 + \w\ 2 — 1 and 2H = C 2 + D 2 = 1 as an additional normalization. The latter means 
that we parametrize geodesies by the natural parameter. Therefore, 

<p = z(A + iB)-w{C + iD), 
ip = z(C + iD)+w(A + iB). 

Let us introduce an auxiliary function p = w/z. Then substituting tp and ip in the Hamiltonian system 
we get the equation for p as 

p= (C + iD)p 2 -2iBp+{C -iD), p(0) = 0. 

The solution is 



(C-iD) sm(sVmP) 

p{s) = 



Vl + B 2 cosisVT+B 2 ) + iB sin(sVl + B 2 ) 
Taking into account that zz = —ww, we get the solution 



?(s) = (cos(s\/l + B 2 ) + i ■ B sin(s\/i + B T 
V VI + B 2 



e- lBs , (5.1) 



and 

wis) = C + lD sm( s ^l + B 2 )e lBs . (5.2) 
Vl + B 2 

If B = we get the solutions with constant horizontal velocity coordinates 

z(s) = coss, w(s) = (2:3(0) + ix4(0)) sins 

from the previous section. 
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Theorem 2. Let A be a point of the vertical line, i. e. A = (cos oj, sin a;, 0, 0), lo £ [0, then there are 

countably many geometrically different geodesies 7„ connecting O = (1,0,0,0) with A. They have the 
following parametric equations 

z n (s) = (cos(s 7 2 )-z-sin( s ™ f ))eVPfe? (5.3) 
= (A 3 (0) + »»4(0)) sin (s 7 - ) e 7?fer, 



V7T 



n£Z \ {0, ±1}, s G [0, s„], where l n = -j^Sn — V^n 2 — cj 2 zs t/ie length of the geodesic j n . 

Proof. Since we use the condition 2_ff = |i| 2 + |w| 2 = 1 we conclude that the geodesies are parametrized 
by the arc length, and the length of a geodesic at the value of the parameter s — ly2 is equal to Z. If 
the point A = (z(s),w(s)) belongs to the vertical line starting at O = (1,0,0,0), then \z(s)\ — 1 and 
|w(s)| = provided that — Bs = lo. It implies 

cos 2 (s\/l + B2 )+ 1 + B 2 si n 2 0vl + B 2 ) = 1, sin(sA/l + B 2 ) = 0, -Bs = u. 
These equations are satisfied when 

s 



= ^v 2 n 2 -uj 2 , Bn = -^===, neZ\{0}. 

We conclude, that for n(s) S Z\{0} there is a constant B n (s) = — g " 2 , such that the corresponding 
geodesic 7„(s), s S [0, s n ], satisfying equation (|5.3p joins the points O and A and the length of the 
geodesic is equal to s„ — \J ir 2 n 2 — lu 2 . □ 

Remark 2. In the formulation of the theorem the words 'geometrically different' mean that due to the 
rotation of the argument of C + iD in w(s), there exist uncountably many geodesies. 

So far we have had a clear picture of trivial geodesies whose velocity has constant coordinates. 
They are essentially unique (up to periodicity). The situation with geodesies joining the point (1, 0, 0, 0) 
with the points of the vertical line A has been described in the preceding theorem. Let us consider the 
general position of points on S 3 . 

Theorem 3. Given an arbitrary point (zi,u>i) € § 3 which neither belongs to the vertical line A nor to 
the horizontal sphere S 2 , there is a finite number of geometrically different geodesies joining the initial 
point (zq,wq) £ § 3 with (zi,wi), zq = 1, wq = 0. 

Proof. Let us denote 

W\ = pe , z\ = re , C + iD = e . 
Then from ([5~Tjl and ([572]) we have that 

r 2 = 1 - - sin 2 (s a/I + B 2 ), and <p = Bs + 9, (5.4) 

where s is the value of the length arc parameter when the point (zi,wi) is reached. We want to exclude 
the parameter s and rewrite the equations (|5.1[) and (|5.2[) in terms of the parameter B and the given 
dates r,p ^ 0,a / 0, and ip. We suppose for the moment that the angles syl + B 2 and sB are from 
the first quadrant. Other cases are treated similarly. Then we have 

z = ( y/l - (1 + B 2 )p 2 + iBp)e l(e ~ v \ 

and 

9 = 9(B) = a + in — arctan — = 

Vi-(i + s 2 V 2 

The first expression in (|5.4[) leads to the value of the length parameter s at (z\,wi) 

s = ; arcsin(p\/ 1 + B 2 ), 

VI + B 2 



8 



Der-Chen Chang, Irina Markina and Alexander Vasil'ev 



and the second to 

tp = 6 H , arcsinfpv 1 + B 2 ). 

VI + B 2 ' 

Substituting 0(B) in the latter equation we obtain 




sin (a -arctan (-5 — ^))v 1 + ^ = p\Jl + B 2 , (5.5) 



as an equation for the parameter B. Observe that tp—9(B) — a — arctan ( ^rrg^ ^ is a bounded function 

and linis^o 7^ 0. Indeed, if the latter limit were vanishing, then the value of given (p would be zero 
and the solution of the problem would be only B — which is the trivial case excluded from the theorem. 
So the left-hand side of equation (|5.5|) is a function of B which is bounded by 1 in absolute value and 
fast oscillating about the point B — 0. Observe, that a = corresponds to the horizontal sphere which 
also was excluded from the theorem. The right-hand side of ()5.5|) is an even function increasing for 
B > 0, see Figure (jTJ) . Therefore, there exists a countable number of non- vanishing different solutions 

{B n } of the equation (|5.5p within the interval \B\ < J — 1 = j^j with a limit point at the origin. 




Figure 1. Solutions to the equation (|5.5 



However, in order do define the parameters B, we need to solve the equations (|5.4|) , (|5.5p . and not 
all B n satisfy all three equations. Let us consider positive B n . We calculate the argument of z as 



-Bri-s + arctan 



= tan (s^BjTl) 



—B n s + arctan 
< -B n s 



B n p 



^\-{l + Bl)p 2 

B n P 



^l-(l + B 2 )p 2 



On the other hand, we have 



1 _ arcsin(p v /l + B 2 X ) p 



syTTBj s 

Observe that due the remark before this theorem, a > and < p < 1. Therefore, we deduce the 
inequality 

l-Vl-p 2 (l + £2) 



a < B n p- 



^/l-p 2 (l + B 2 n ) 
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Vl-^(l +jB 2) 

B„p>a v ==. (5.6) 

1- Vl-p 2 (l + ^2) 

The right-hand side of the inequality (|5.6p decreases with respect to -B„ > 0. 

Set e = ^j^- If e < p 2 (l + -B 2 ) < 1, then immediately we have the inequality B 2 > - 1) > 0. 
If < p 2 (l + B 2 ) < e, then the inequality (|5.6|) implies that 



B„ > a— . = a — > 0. 

p(i-vr^i) p(V2-v/r~^) 

Finally, we obtain 



This proves that all positive solutions to the equation (|5.5p must belong to the interval (6(£i, p), y — 1), 
hence there are only finite number of such B n . The same arguments are applied for negative values of 
B n . 

Let us discuss the limiting cases. If p — * then the endpoint aproaches the vertical line. In this 
case the graph of the right hand side function in (|5.5[) approaches the horizontal axis and the range of 
\B\ increases. If a — > then the end point aproaches the horizontal sphere S 2 , the number of geodesies 
is finite for any value a ^ 0, but decreases. □ 

This theorem reveals similarity of sub-Riemannian geodesies on the sphere with those for the 
Heisenberg group. The number of geodesies joining the origin with a point neither from the vertical 
axis nor from the horizontal plane is finite and approaching the vertical line becomes infinite. 

Hyperspherical coordinates 

Let us use the hyperspherical coordinates to find geodesies with non-constant velocity coordinates. 

x\ + ix 2 = e^ 1 cosr], (5.7) 
x 3 +ix 4 = e^ 2 sinry, ?/e[0,7r/2), £i, £ 2 S [-7T, tt) 

The horizontal coordinates are written as 

sin 

a = r)cos(£i - £ 2 ) + (6 + £ 2 )sin(£i - £ 2 X 



b = -r)sin(£i - £ 2 ) + (£i + £ 2 ) cos(£i - £ 2 ) 
c = £i cos 2 77 — ^2 sin 2 77. 



2 ' 
sin 2 77 



The horizontality condition in hyperspherical coordinates becomes 

f 1 cos 2 77 — £2 sin 2 7] = 0. 

The horizontal sphere (|4.4| is obtained from the parametrization (|5.7p . if we set £1 = 0, £ 2 = ip, 
i] = s. We get 

a 2 + b 2 = 1 = f] 2 =>• a = cost/;, — sin-0. 

The vertical line is obtained from the parametrization (|5.7p setting 7/ = 0, £1 = s. 
Writing the vector fields N, Z, X, Y in the hyperspherical coordinates we get 

N = — 2 cotan 277(9,, , Z = — d^ 2 , 
X = sin(£i - £ 2 ) tan + sin(£i - £ 2 ) cotan 77% + 2cos(£i - &)d v , 
Y = cos(^i - £ 2 ) tan 77% + cos(£i - £ 2 ) cotan 77% - 2sin(£i - £ 2 )(V 
In this parametrization the similarity with the Heisenberg group can be shown. The commutator of 
two horizontal vector fields X, Y gives the constant vector field Z which is orthogonal to the horizontal 
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vector fields at each point of the manifold. In hyperspherical coordinates it is easy to see that the form 
to = cos 2 — sin 2 77^61 that defines the horizontal distribution is contact because 

oj A duj = sin(2?7) dr] A d£i A d& = 2dV, 

where dV is the volume form. The sub-Laplacian is defined as 

i(X 2 + Y 2 ) = i(tan 2 r)dl + cotan 2 ryc»| 2 + 4d 2 + 2%%). 
The Hamiltonian becomes 

H(£i,&,r),i>i,i>2,0) = 2 ( tan2 *fil>i + cotan2 #1 + 46,2 + ^2 ) , 
and the corresponding Hamiltonian system is given as 

Si = 77— = tan 77 + 

; 9H 2 

6 = 77— = ^2 cotan z 77 + t/>i 

<96 
#6 

• 9 tan 77 , o cotan 77 

or/ cos 77 sin 77 

Let us solve this Hamiltonian system for the following initial data: 77(0) = 0, 6(0) = 0, 6(0) = Oj 

Mo) = M Mo) = M o(o) = = e . 

We see that t^i and ip2 are constant. The horizontality condition at (0, 0, 0) gives 60) = and the 
first equation of Hamiltonian system implies that ip2 = 0. If ipi = 0, then 6 = 6 = 0, 77 = 46*o and we 
get the variety of trivial geodesies (|4.4[) up to the reparametrization s 1— ► A6$s. To find other geodesies 
we suppose that / and 77(C)) > 0. This condition ensures us that the trajectory starting at the 
point (0, 0, 0, ) remains in the domain of parametrizaition locally in time. From the third and from the 
last equations of the Hamiltonian system we have 

A ,2 Sm? 7 ■ T A ,2 Sm? 7 7 

f] = —Aip 1 — => T]dr] — ~Atp 1 — d?7 

COS 77 COS 77 



t) 2 =C 



9 1 

cos z ?y 



We observe that C = t) 2 (0) + Aip\ > 0. 

Continue to solve the Hamiltonian system finding 77(5) 

cos 77 dr] 



Denote by sin 77 = p. Then, 

Integrating (|5.8[) from to s we get 

s + K-- 



\J C cos 2 77 — 
dp 



ds. 



rcsin (\/^5f sinr ' (s) ) 



= ds. (5.8) 
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where K is found setting s = as K = —i — arcsinO = 0. We calculate 



sin 2 n(s) = ° ^ sin 2 (VCs). (5.9) 
From the Hamiltonian system we find 

600=Vi«, (5.10) 

and 

• _ sin 2 77(a) _ sin 2 (VCs) _ C 

$1 ~ % - sin 2 V (s) ^ 1 ffl + sin 2 (v^ S ) ' ° " C^ 2 ' 

It gives 



£i(s) = — t/»is H — py-r arctan '^l 1 tan(\/Cs) . (5.11 



2|V 



2|^i| L 

Let us suppose for the moment that the geodesies are parametrized on the interval [0,1]. If the 
initial point and the finite point are on the vertical line: 77(0) = 77(1) = 0, then 

= sin 2 77(1) = ^-Msin 2 (VC) C = ttV. 

O 

Since the value of £2 on the vertical line is arbitrary, the values of C and £i(l) give us the value of £2(1) ■ 
Setting C = 7r 2 n 2 in the equation for £i(l), we find ipi = — £i(l). Then 

6(i) = -£i(l). 

The finite point on the vertical line corresponds to the value of £i(l),?7(l) = 0, and £2 = — 

We also note that the square of the velocity \v\ 2 = rj 2 (s) + (£i(s) + ^2(s)) 2 s '" ^ s constant 

along geodesies. Applying the initial condition 77(0) = 0, we get 

| V | 2 = 7) 2 (0) = C-4^ 2 . 

In the case when a geodesic ends at the vertical line at £i(l), its lengths is expressed as 



V2l n = sjC-A^l = ^/W-4£(l). 

We see that this result coincides with the result given by Theorem [2] and we state it as follows. 

Theorem 4. Let A be a point of the vertical line, i. e. £1 is given and 77 = 0. There are countably many 
geodesies connecting O — (0,0,0) and A, given parametrically as 

£i( s ) — £i s + « arctan ( — - tan(7rns)V 

600 = -£is, 

shirks) = — sin(7rns), 



where -j^s n = -^V 7r2n ' 2 — ^£1 * s ^ e ^ en gth of the geodesic j n , n S 



6. Hopf fibration 

There is a close relation between the sub-Riemannian sphere S 3 and the Hopf fibration. Let § 2 and § 3 
be unit 2-dimensional and 3-dimensional sphere respectively. We remind that the Hopf fibration is a 
principal circle bundle over two-sphere given by the map h : § 3 — > § 2 : 

h(xi,X2,xa,X4,) = {{x\ + x\) - {x\ + x\), 2{xiXi + x 2 x 3 ), 2(x 2 x i - xix 3 )). 

Another way to define the Hopf fibration is to write 

h(q) = qiq* G § 2 , g £ § 3 , i = (0,1, 0,0). 

The fiber passing through the unity of the group (1, 0, 0, 0) has equation (cos#, sin#, 0, 0), which as we 
see, coincides with the equation of the vertical line at this point. The sphere S 2 represents the horizontal 
"plane" sweep out by the geodesies with constant horizontal coordinates. 
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Definition 1. Let Q — > M be a principle G-bundle with the horizontal distribution T> on Q. A sub- 
Riemannian metric on Q that has distribution 2? and it is invariant under the action of G is called a 
metric of bundle type. 

In our situation S 3 — ► § 2 is a principle S^-bundle given by the Hopf map. The sub-Ricmannian 
metric on the distribution T> = span{A, Y} was defined as the restriction of the euclidean metric (•, •) 
from R and we used the same notation (•, •) for sub-Riemannian metric. 

Proposition 3. The sub-Riemannian metric (•, •) on § 3 is a metric of bundle type. 

Proof. The action of the group 8 1 on q= (xi,X2,X3,%4) £ § 3 can be written as qoe lt , e lt — (cost + 
isint) £ S 1 , t £ [0,2%), where o is the quaternion multiplication. If we write q = (e** 1 cos 77, e^ 2 sin rj), 
q = qo e lt , then q — (e l ^ 1+t ) cos 77, e 1 ^ 2- '-* sin rj) . In order to show that the metric (•, •) is ofbundle type, 
we have to prove that the metric is invariant under the action of the group § . The metric (•, •) at any 
q is given by the matrix 

"1 " 

cos 2 r\ 
sin 2 rj 

and it is easy to see that it is invariant under the action q — q o e lt . □ 

We can formulate the results of Theorems [2] and [4] as an isoholonomic problem. First let us give 
some definitions. Let c : [0, 1] — > § 2 be a curve in § 2 For a given point x of the fiber Q c (o) at c(0), 
let 7 be a horizontal lift of c that starts at q (it means that the projection of 7 under the Hopf map 
h : § 3 — > § 2 coincides with c). The map <I>(c) : Q c (p) — > Qc(i) sending x — 7(0) to the endpoint 7(1) 
of the horizontal lift is called parallel transport along c. The action of S 1 takes horizontal curves to 
horizontal curves, so any two horizontal curves 71 and 72 of c are related by ji = 72 <? for some g £ S 1 . 
It follows that the action of S 1 commutes with $(0), that is, $>(c)(xg) = ($(c)(x))g. 

If c is a closed loop, parallel transport $(0) maps the fiber Q c (o) onto itself. Fix a point xq £ Q c (o)- 
Since S 1 acts transitively on Q c (q) we have <p(c)(xo) = 2:0^ for some I £ 8 . If we choose another point 
Vo = x g £ Qc(o) , we get 

$(c)(y ) = (*(c)(a ))(ff) = x lg = y^g^lg). 

The curve c therefore, determines a conjugacy class in S , called the holonomy class of c. The element 
/ £ S 1 for which <i>(c)(xo) = xqI is called the representative holonomy of c with respect to xq. The set 
of all such I £ S 1 for c running over all closed loops with c(0) = c(l) = ft.(xo) is a subgroup of S 1 called 
the holonomy group of the distribution T> &t Xq. 

Let us fix a representative holonomy I £ S 1 and a point xo € § 3 , or equivalently the initial and 
the finite points Xq and X\ = XqI. The set of horizontal curves that join xq to xi is in one-to-one 
correspondence with the set of all closed loops on S 2 , based at h(xo), whose holonomy with respect 
to xo is I. Recall that the Riemannian length of a loop on § 2 equals the sub-Riemannian length of its 
horizontal lift. Thus, the sub-Riemannian geodesic problem for geodesies with endpoints at the same 
fiber is equivalent to the following isoholonomic problem: Among all loops with a given holonomy, find 
the shortest. 

If we take xq = (1,0,0,0), I = e lu} , x\ = x$l = (cosw, sinw, 0, 0), then Theorem [2] says 
if ui £ [0, tt) then the shortest loop has length si = \/ir 2 n 2 — uj 2 , 
if uj £ [tt, 2ir) then the shortest loop has length S2 = V 47r 2 n 2 — lo 2 . 
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